Multiplicity fluctuations due to the temperature fluctuations in high-energy nuclear 

collisions 
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We investigate the multiplicity fluctuations observed in high-energy nuclear collisions attributing 
them to intrinsic fluctuations of temperature of the hadronizing system formed in such processes. 
To account for these fluctuations we replace the usual Boltzmann-Gibbs (BG) statistics by the non- 
extensive Tsallis statistics characterized by the nonextensivity parameter q, with \q — 1| being a 
direct measure of fluctuations. In the limit of vanishing fluctuations, q — » 1 and Tsallis statistics 
converges to the usual BG. We evaluate the nonextensivity parameter q and its dependence on 
the hadronizing system size from the experimentally observed collision centrality dependence of the 
mean multiplicity, (TV), and its variance, Var(N). We attribute the observed system size dependence 
of q to the finiteness of the hadronizing source with q = 1 corresponding to an infinite, thermalized 
source with a fixed temperature, and with q > 1 (which is observed) corresponding to a finite source 
in which both the temperature and energy fluctuate. 

PACS numbers: 25.75.Ag, 24.60.Ky, 24.10.Pa, 05.90.+m 



I. INTRODUCTION 

With the large number of particles produced in heavy 
ion collisions at CERN SPS and BNL RHIC it is possi- 
ble to study fluctuations in different physical observables 
on the event-by-event basis p]]. These fluctuations are 
potentially a very important source of information on 
the thermodynamical properties of strongly interacting 
systems formed in such collisions, like, for example, its 
specific heat |2| (connected with fluctuations observed in 
particle multiplicities [H, HJ], in transverse momenta @ 
and in other global observables), its chemical potential 
or matter compressibility @. 

Fluctuations of multiplicity observed in heavy ion col- 
lisions exhibit spectacular and unexpected features as 
functions of the number of participants. Recent results 
on centrality dependence in Pb + Pb collisions at 158A 
GeV obtained by NA49 and WA98 Q experiments 
indicate that the scaled variance of the multiplicity dis- 
tribution, Var(N)/ (N), increases when proceeding from 
the central towards peripheral collisions, i.e., when the 
number of participants decreases. Such behavior is con- 
firmed by a comprehensive survey of multiplicity fluctu- 
ations of charged hadrons provided by PHENIX experi- 
ment 0- 

During the last decade the statistical models of strong 
interactions were constantly used as an important tool to 
study fluctuation pattern observed in high energy nuclear 
collisions experiments, which are quantified by scaled 
variance mentioned above (see, for example, [a 0, [l(| 
and references therein). However, to minimize the effect 
of the possible participant number fluctuations, the anal- 
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ysis of particle number fluctuations has been restricted 
to the most central A + A collisions only. 

However, up to now, none of the models aimed to de- 
scribe all essential features of multiparticle production 
processes (like multiplicities and distributions of particles 
in phase space and their composition) describe also the 
dependence of Var(N) / (N) on Np, the number of nucle- 
ons from projectile nucleus participating in the collision - 
projectile participant, observed experimentally Q. They 
lead to multiplicity distributions of the (approximately) 
Poissonian form, independent of centrality, which highly 
underestimate the observed multiplicity fluctuations in 
noncentral collisions and thus are unable to reproduce, 
even qualitatively, the observed centrality dependence of 
the scaled variance. This remark applies not only to 
the Monte Carlo models like H1JING [n|, HSD [l| or 
UrQMD [l3[, which are based on string excitation and 
decay, but also to any statistical model that does not 
assume correlations among secondary particles (cf. dis- 
cussion in Sec. IIII Cp . On the other hand, these results 
can be described by some specialized models addressing 
fluctuations directly, like the percolation model the 
model assuming inter-particle correlations caused by the 
combination of strong and electromagnetic interactions 
[IH or the transparency, mixing and reflection model [l6| . 
Fluctuations in these models reflect some dynamical fea- 
tures of the production process, (specific for the model 
considered). 

In this paper we shall address the problem of multiplic- 
ity fluctuations without resorting to any specific dynami- 
cal picture but, instead, by attributing them to some non- 
statistical, intrinsic fluctuations existing in a hadronizing 
system produced in high energy heavy ion collisions. To 
account for such fluctuations we shall use a special ver- 
sion of statistical model based on nonextensive Tsallis 
statistics [TtJ in which fluctuations of the temperature are 
known to be directly connected with the nonextensivity 
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parameter q [HHH, namely q = 1 + Var(l/T)/(I/T) 2 . 
In what follows they will be assumed to be the true (if 
not the only) origin of the fluctuations observed in the 
experimental data. The resulting distributions are then 
power-like, exp(-E/T) =>■ exp q (-E/T) = [1 - (1 - 
q)E/T] 1 /( 1 ~'fi and in the limit of q —> 1 they go smoothly 
to the usual Boltzmann distributions (in what follows we 
shall address only the q > 1 case) 20] . It is important for 
our further discussion that such non-exponential distri- 
butions of energy result in a non-Poissonian multiplicity 
distributions of the produced secondaries [25| . 

However, to incorporate the new features of data re- 
ported by [20|, we have to extend the notion of fluctu- 
ating temperature replacing it by some g-dependent ef- 
fective temperature, T e tf, which accounts not only for 
the intrinsic fluctuations in the hadronizing source (as in 
[HI]) but also for effects of the possible energy transfer 
taking place between the hadronizing source and its sur- 
roundings {27j. This will be done in Section [TTJ Section 
IIIII contains our results: the universal participant depen- 
dence (i.e., scaling in the variable / = Np/A) presented 
in Section fill Al its explanation by the observation that 
in reality the hadronizing source is always of finite size 
presented in Section UlI Bl and the system size dependence 
of multiplicity and multiplicity fluctuations presented in 
Sections IIII CI Section IIVI summarizes our work. Some 
details are presented in Appendices and [Bj 




II. EFFECTIVE TEMPERATURE 

In the proposed approach we replace the standard 
Boltzmann-Gibbs exponential distribution, 

g(E)=Cexp(-E/T), (1) 

by the Tsallis distribution (^-exponential) defined as 



FIG. 1: (Color on line) Dependence of the effective tempera- 
ture T e f f (in GeV) on the parameter q for the production in 
different reactions of negative pions (o) and antiprotons (b), 
all data points are from [2y| . The solid lines show linear fits to 
the obtained results: for upper panel T e // = 0.22 — 1.25(g— 1), 
for lower panel T e // = 0.36 — 3A(q — 1). 
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with c p , p and D being, respectively, the specific heat 
under constant pressure, density and the strength of the 
temperature fluctuations (cf. Appendix lAl for details). 
Effective temperature T e ff occurs when we experience 
both the fluctuations of the temperature T (around the 
value T ) and some energy transfer taking place between 
the source and the surroundings given by <f> [28[. Notice 
that in Eq. ([!]) energy transfer affects T e ft only in the 



presence of fluctuations, i.e., for q ^ 1. The predicted 
q dependence of T e f f is indeed observed experimentally, 
cf. Fig. Q] Namely, in [26| the transverse momentum 
spectra of pions and antiprotons produced in the inter- 
actions of p+p, d+Au and Au+Au at ^snn = 200 GeV 
at RHIC-BNL experiments [30| were analyzed using a 
nonextensive approach in which the slope of px distribu- 
tion determines the effective temperature T e ff. Its shape 
gives the parameter of nonextensivity q (3lj | . They evalu- 
ated, among other things, the nonextensivity parameter 
q and the effective temperature T e f f for a different num- 
ber of participants, Np. From them we can deduce the 
dependence of T e f / on the parameter q which is shown 
in Fig. [TJ Notice that in all cases we find that (j> < 
and that T e ff seems to depend linearly on (q— 1). Nega- 
tive values of <j) mean that the energy is transferred from 
the interaction region to the surroundings (i.e., to the 
spectators of noninteracting nucleons) [33j |. 
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III. DIFFERENT FACETS OF MULTIPLICITY 
FLUCTUATIONS 

A. The universal participant dependence 

We start a discussion of multiplicity fluctuations by 
recollecting the result of [25[ saying that if N particles are 
distributed in energy according to the ./V-particle Tsal- 
lis distribution described by the nonextensive parame- 
ter q then their multiplicity distribution has to be of the 
Negative-Binomial type with fc _1 = q — 1 (cf. Appendix 
|B]for details). It means therefore that we can expect that 



Var(N) _ -, 

(N) 



(5) 



(and, what is important, that q — 1 determined this way 
does not depend on the acceptance of the detector [3l|). 
On the other hand, if U is the accessible energy and gT e f / 
is the mean energy per particle detected in the acceptance 
region (with g being a parameter and T e f / effective tem- 
perature defined in Eq. Q ), than 



(N) 



(U) 

9 T eff ' 



Using Eq. (j4j one can shown that 



(N) - nt) Np 
(N) 



= c(q-l), 



(6) 



(7) 



where no is the multiplicity in the single nucleon-nucleon 
collision measured in the region of acceptance, it is de- 
fined by the constraint that (U) — n^gT^Np, whereas 
c = —(f)/ (Dc p pTo) is a constant (notice that because in 
cases we are interested here (f> < 0, cf. Fig. [IJ c is pos- 
itive). Comparing now Eqs. © and (HJ) one can expect 
that 



Var(N) 



l + c((N) -n N P ) 



Notice that Eqs. ([7]) and ([5]) do not depend on parameter 
g anymore; its role in present work is to fit, if necessary, 
(N) in Eq. to experimental data only. 

This means that using the notion of T e ff one should 
observe that Var(N) / (N) and (N) are mutually con- 
nected. Fig. [2] shows that this is indeed the case. If 
T e ff depends on q (i.e., for c ^ 0) than the dependence 
Var(N)/ (N) on the number of participants Np is con- 
nected to the dependence of (N) on Np. Notice that if 
(N) is linear in Np then Var(N)/(N) is constant. In 
particular, for (N) = n N P we have Var(N)/(N) = 1. 
Therefore the experimental fact that Var(N) / (N) de- 
creases with increasing Np indicates nonlinear depen- 
dence of (N) on the number of participants Np. 

It should be stressed that data taken by the NA49 
experiment @ show that Var(N)/ (N) changes rather 
strongly with Np. For peripheral collisions (i.e., for small 
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FIG. 2: (Color online) (a) Comparison of Var (N)/{N) versus 
Np (squares) with 1 + c({N) — noNp) versus Np (circles) 
for no = 0.642 and c = 4.1 37]- Data are for negatively 
charged particles from Pb+Pb collisions as collected by NA49 
experiment [■( . (b) The same but translated to dependence of 
q— 1 versus Np . Squares were obtained from the Var(N)/ {N) 
versus Np dependence whereas circles from the (N) versus Np 
dependence. The values of no and c are the same as before. 
The solid line shows the dependence (|16p for the nuclear mass 
number A = 207 and parameter a = 0.98. 



Np) one observes marked deviation of Var(N) / (N) from 
unity. At the same time deviation of (N) from linearity in 
Np is very weak. Our result based on the concept of effec- 
tive temperature, T e //, and showing that Var(N)/ (N) 
and (iV) are connected is therefore by no means trivial. 

We close this Section by reminding that all previous 
attempts in this field were addressing only the behav- 
ior of Var(N)/ (N) versus Np, leaving aside the possible 
nonlinear dependence of ( N) on Np . Fig. O shows that 
introducing T e f f results in a kind of scaling, namely the 
same value of q — 1 describes dependencies of both vari- 
ance and mean multiplicity on Np, following Eq. ([5]), 
which - we repeat - has its origin in the ^-dependence of 

TeffM- 
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B. Finite hadronizing source size and temperature 
fluctuations 



We shall now derive the q—1 dependence on Np seen 
in Fig. [5] Let us first observe that 0, |39| 



Var(T) 



c\. 



(9) 



i.e., the parameter q can be regarded as connected (via 
fluctuations of temperature) to the heat capacity under 
constant volume, Cy. For a system with finite size re- 
maining in contact with a heat bath one has, following 
Lindhard's approach [4fj], that 



Var(U) + C v Var(T) = (T) 2 C V , 



(10) 



This is a kind of uncertainty relation (in the sense that 
it expresses the truth that in the case of conjugate vari- 
ables one standard deviation in some measurement can 
only become small at the expense of the increase of some 
other standard deviation [41] ) . Relation (fT0|) is supposed 
to be valid all the way from the canonical ensemble, 
where Var(T) = and Var(U) = (T) 2 C V , to the micro- 
canonical ensemble, for which Var(T) = (T) 2 /Cv and 
Var(U) — 0. Eq. (jTUJ) expresses both the complemen- 
tarity between the temperature and energy, and between 
the canonical and microcanonical description of the sys- 
tem pf. 

To obtain realistic (intermediate) distributions, start 
from a system at a fixed temperature T. The standard 
deviation of its energy is 

Var(U) = <T) 2 ^ = (T) 2 C V . (11) 
Inverting the canonical distribution Qt{U) one can obtain 



gu(T) = -T 2 



d_ 

dT 



g T {U')dU' 



(12) 



and interpret it as a probability distribution of the tem- 
perature in the system (cf., Eq. (|A10|) ). The standard 
deviation of this distribution then yields [44[ 



Var(T) = 



C v 



(13) 



Because for a canonically distributed system the energy 
variance is Var(U) — (T) 2 Cv and because for an iso- 
lated system Var(U) = 0, for the intermediate case the 
variance (expressing energy fluctuations in the system) 
can be assumed to be equal to 



Var(U) = (T) 2 C V ^ £€(0,1), 



(14) 



where the parameter £ depends on the size of the 
hadronizing source. Inserting this in Eq. (JTUJ) one gets 
that q depends on £ in the following way: 



Var(T) 1 - £ 



(15) 



Assuming now that the size of the thermal system pro- 
duced in heavy ion collisions is proportional to the num- 
ber of participating nucleons Np, i.e., that £ ~ / = 
Np/A, and taking into account that Cy = aNp, we ob- 
tain that 



1 



1 



aN, 



■(I-/)- 



(16) 



This is the relation which nicely fits data, cf., Fig. [2] 

To recapitulate: we know that if U = const and 
T = const then the multiplicity distribution P(N) is 
Poissonian (cf. Appendix IB 1[| and we also know how 
fluctuations of T change P(N) from the Poissonian form 
to the NBD one (cf . Appendix IB 2[) . We want now to see 
how big are fluctuations of T in our hadronizing systems 
formed in a collision process. It turns out that for N — 
const we have either Eq. (fTTj) or Eq. (fl"3|) depending on 
whether T = const or U = const. We claim then that 
we can assume validity of Eq. (fT0|) not only in the above 
limiting cases but also in the general case when both the 
energy U and the temperature T fluctuate at the same 
time, i.e., if fluctuations of the energy U are given by Eq. 
(fTl"]) then fluctuations of the temperature T are given by 
Eq. Q15p. Finally, knowing how big are fluctuations of T 
we can deduce fluctuations of the multiplicity N. This 
will be performed in what follows. 



C. System size dependence of mean multiplicity 
and multiplicity fluctuations 

We shall now discuss, respectively, the system size de- 
pendence of the mean multiplicity and multiplicity fluc- 
tuations. In our approach the system size enters through 
Cy and £. We shall keep £ = Np/A and connect Cv 
with the measured quantities considering two natural as- 
sumptions: either 



or 



C v = aNp 



C v - a'(N). 



(17) 



(18) 



When used together with Eqs.([7|) and (fT5]) the second 
possibility leads to a very simple scaling relation: 

W-§>)|^ = £(l-^) (19) 
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whereas the first one results in a slightly more involved 
formula 



(20) 



In both cases (N)\n p =a = tiqA is multiplicity extrapo- 
lated to Np — A. As seen in Fig. 02 where compari- 
son with experimental data [|[ is presented, one observes 
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different dependencies for Pb + Pb collision and lighter 
nuclei for which in scmiccntral collisions such an effect is 
practically not observed. In addition, for peripheral colli- 
sions (Np < 0.15^4 ) one observes the deviation from the 
expected dependence (fT9|) . Notice that deviation from 
the linear fit for Pb+Pb collisions concern only the 5 most 
peripheral points and the observed discrepancy means 
that the measured experimentally mean multiplicity is 
less than 1 particle higher than the expected value. The 
agreement with prediction given by Eq. (|20p seems to be 
better, what means that the assumption Cy = aNp is 
more realistic. 



following simple scaling relation: 
Var(N) , ., , 1 



1 - 



N P 



(21) 



i.e., in the multiplicity fluctuations dependent on the 
fraction of nucleons participating in the collision / = 
Np/A. Such scaling was found recently by the NA49 col- 
laboration Q when comparing collisions of C+C, Si + Si 
and Pb + Pb, cf. Fig. gl If, instead, we use (JTTJ) then, 
taking into account dependence of (N) on Np given by 
Eq. (JT]), one obtains that, more exactly, 
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FIG. 3: (Color online) Mean multiplicities compared with 
predictions given by Eq.([20} (a) and Eq. (fT9)l (b). We adopt 
n = 0.559, 0.575, 0.657 and 0.642 for, respectively, p + p, 
C + C, Si + Si and Pb + Pb collisions. Linear fit show the 
predictions of Eq. (I19[l with parameter c/a = 2.75 and Eq. 
(|20[) with parameter c/a = 1.8. Only statistical errors are 
indicated. 



Var(N) 



l + (N)(q- !) = ! + ■ 



"o (i - Y 



c ( A 
A\N P 



1 



(22) 



We observe thus a weak dependence on the mass num- 
ber A of colliding nuclei. For a small number of par- 
ticipants, Np « A, we observe an additional increase 
of relative variance, Var(N)/ (N) , in comparison with 
prediction (f2"Tj) . Notice that for <fi > one expects 
just the opposite trend, i.e., the non-monotonic behavior 
of Var(N)/ (N) with increasing number of participants, 
Np. The question about the sign of the heat source term 
(f> seems to still be open due to the lack of data in the 
region of small number of participants. 

For central collisions (i.e., for / = Np/A ~ 
1) multiplicity distributions are sub-Poissonian and 
[Var(N)/(N)]\ Np= A — loq < I. For example, within 
a statistical model with fixed volume loq varies in the 
range 0.5 1.0 [45]. In particular (see [1, [l(| and refer- 
ences therein): (i) The global conservation laws imposed 
on each microscopic state of the statistical system lead to 
suppression of the particle number fluctuations. The final 
state scaled variance behavior in the canonical ensemble 
is characterized by ujq = 0.8 whereas in the microcanoni- 
cal ensemble by ujq = 0.3. (ii) In the primordial values of 
the scaled variance (i.e., calculated before decays of reso- 
nances) one can also observe the effect of quantum statis- 
tics. It turns out that proper inclusion of Bose statistics 
leads to the additional enhancement of particle number 
fluctuations of the order of ujq = 1.05 + 1.06, which is 
quite small at the chemical freeze-out. (Hi) The particle 
number fluctuations can be also enhanced by including 
explicite production and decay of resonances (either in 
the grand canonical ensemble of in canonical ensemble). 
They lead to ujq = 1.1. The actual value of parameter ujq 
depends on the energy of collision and on the acceptance 
in which the multiplicity fluctuations are observed and 
varies in the range 0.8 — 1.0 [46J. For this reason in Fig. 
2] we compare experimental data with the formula 



Var(N) 



1 



N P 



(23) 



In what concerns multiplicity fluctuations, assumption 
(fT5)l results (when combining Eqs. (fTo| and ([5])) in the 



Finally, we note that the transverse-momentum fluc- 
tuations measured in nuclear collisions at 158 A GeV [j| 
and quantified by the measure $ (pt) [47| show a similar 
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centrality dependence, namely, as seen in Fig. [5l 
$ Or) = $n p =a + b(l- -fj . 



(24) 



This behavior of the transverse-momentum fluctuations 
as a function of collision centrality was related in super- 



position model to the centrality dependence of the mul- 
tiplicity fluctuations [49J using the correlation of average 
transverse momentum and multiplicity observed in colli- 
sions @. To derive Eq. (j2~4")l notice that in general $ is 
governed by the multiplicity fluctuations in the following 
way [5(ij |: 



J 



\ 



Var (p T ) + 2(p T ) 



a Var(N) 
(N) 



1 - 



Varipr) (N) 
P \ (pr) 2 Var(N) + 



- y/Var (p T ), 



(25) 



where p is the correlation coefficient between N and 
'Ylii PiT- As in [49j| we can see that multiplicity fluctua- 
tions determine the behavior of <& (pt)- In the first order 
approximation and taking into account multiplicity fluc- 
tuations given by Eq. ([23| we can write the transverse 
momentum fluctuations measure as 



$ (p T ) ~ y/Var (p T ) 



u>o(l - p) 



(Pi 



+ (l-p)y/Var (p T ) 



Var (pt 
(Pt) 2 1 



1 

2 P 



+ 



Var (pt) a' 



1 - 



N P 

T 



,(26) 



which has the form of Eq. ([24]) in what concerns de- 
pendence on the variable Np/A. The slope parameter b 
corresponds to the correlation coefficient p ~ 0.99 (for 
the previously estimated value of a' and for the trans- 
verse momentum fluctuations Var(px)/ (pt) 2 — 0.43, as 
observed experimentally [5l|). Once again we observe 
scaling behavior in the variable / = Np/A. 



IV. SUMMARY 

We have investigated some features of the multiplic- 
ity fluctuations observed in recent high-energy nuclear 
collisions by attributing them to the intrinsic, nonsta- 
tistical, fluctuations of the temperature of the hadroniz- 
ing system produced in such collisions. To this end we 
used the Tsallis statistics approach in which such fluctu- 
ations are accounted for by the nonextensivity parame- 
ter q (more precisely by q — 1). When fluctuations are 
vanishing q — > 1, Tsallis statistics becomes the usual 
Boltzmann-Gibbs (BG) one and the power-like Tsallis 
g-exponents become the usual exponential distributions. 
We were considering a generalized g-exponential ensem- 
ble describing system of Np participating nucleons (as- 
sumed to be proportional to the size of hadronizing sys- 
tem). It was found that in this case one can associate the 
nonextensivity parameter q with the number of partici- 
pants Np (cf. Eq. (Til)]0 . It was found that data require 
a generalization of the usual notion of fluctuating tem- 
perature by adding effects of the energy transfer between 
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FIG. 4: (Color online) The scaled variance of the multiplicity 
distribution of negatively charged particles produced in p+p, 
semicentral C+C, semicentral Si+Si, and Pb+Pb collisions 
as function of the fraction of nucleons participating in the 
collision, Np /A. The experimental data [2| are compared with 
prediction (|23p with parameters ujq — 0.79 and a' — 1.0 (solid 
line) and with prediction (|22[) with parameters luq = 0.87 and 
a = 0.83 (dashed line, the other parameters are the same as 
in Fig. 0with n = 0.642 and c = 4.1 for A = 207). 



the hadronizing source and the surroundings composed 
with nucleons not participating directly in the reaction. 
This resulted in the introduction of a g-dependent effec- 
tive temperature T e ft. We have also allowed for fluctu- 
ations of the full accessible energy U (cf., Section lill A[) . 
without these fluctuations Var(N) / (N) would be simply 
constant, independent of Np/A. 

The simplest possible explanation of the observed ef- 
fect was the q — 1 dependence on the number of partic- 
ipants Np provided by the Eq. p6)l . We assumed that 
parameter connected with the size of the collision region, 
£, is approximately given by £ = Np/A. The modifica- 
tions of this assumption are possible and lead to better 
agreement with experimental data. Among others, the 
presence of nonlinear terms in the £ variable in Eq. (|15|) 
would approve consistency with data (or, for example, 
assuming that £ = ((Np/A) with ((Np/A) more compli- 
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FIG. 5: (Color online) Transverse momentum fluctuations of 
negative particles defined by the measure "3? (4f| as a func- 
tion of number of the fraction of nucleons participating in the 
collision, Np/A. The triangles correspond to p + p collisions, 
asterisks to C + C, squares to Si + Si, and circles to Pb + Pb. 
The experimental data are taken from [f| and are compared 
with prediction of Eg . ([24)l for parameters 3>jv P =A = 0.48 MeV 
and b = 4.8 MeV. 



ity parameter q and its dependence on N p obtained from 
comparison with experimental data does nor depend on 
experimental acceptance. 
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APPENDIX A: TEMPERATURE 
FLUCTUATIONS 

We shall collect here the main points concerning the 
idea of temperature fluctuations Qj3, [27| . 

Suppose that we have a nonhomogeneous thermody- 
namic system in different regions of which there is differ- 
ent temperature T, which fluctuates around some mean 
temperature To. As result of these fluctuations, the ac- 
tual temperature T equals 



cated as used here where ((N P /A) — Np/A). Neverthe- 
less, we restricted our discussion to the simplest possible 
approximation to demonstrate the connection between 
temperature fluctuations and observables from relativis- 
tic ion collisions in a more transparent way. 

Let us close with some remarks concerning limitations 
of our analysis. Firstly, let us first remind that multiplic- 
ity fluctuations considered in this study were observed 
for a fixed number of the projectile participants. Al- 
though, in the collisions of identical nuclei, the average 
numbers of the participants from the projectile equal ap- 
proximately that from the target, the actual number of 
target participants fluctuate whereas the number of the 
projectile participants is kept fixed. However, one should 
notice that the multiplicity fluctuations were measured 
in the forward hemisphere (1.1 < y c . m , < 2.6) in which 
the influence from the target participants is marginal. 
The models which produce, independently of central- 
ity, approximately Poissonian multiplicity distributions 
highly underestimate the observed multiplicity fluctua- 
tions in noncentral collisions. They are unable to re- 
produce, even qualitatively, the centrality dependence of 
the scaled variance (cf., for example, Q and references 
therein) . The highest scaled variance is obtained from the 
hadron-string dynamics [l2j and ultrarelativistic molec- 
ular dynamics [13|. However, both approaches (in their 
standard versions) show a flat scaled variance, u> ~ 1.2, 
and exhibit almost no dependence on N p [54j. Finally, 
note that in comparison with experimental data (which 
have finite acceptance, equal to p = 0.16 in the case of 
NA49 data Q ) we have used multiplicity N given by this 
acceptance. Also all parameters used for description of 
these data (like no, a', ojq, &n —a an d b) are acceptance 
dependent, However, as shown in [36| . the nonextensiv- 



T = To - £(t)T, 



(Al) 



where £(t) describes the actual (not specified) stochastic 
process causing these fluctuations. 

There must take place an exchange of energy (heat) 
between the regions mentioned above, in particular be- 
tween any selected region and the rest of the system. This 
exchange eventually leads to equilibration of the temper- 
ature of the whole system. The corresponding process of 
heat conductance is described by the following equation 



c p p- 



dT 
dt 



7(1" — T) = 4>, 



(A2) 



where c p , p and 7 are, respectively, the specific heat un- 
der constant pressure, the density and the coefficient of 
external conductance. The heat source term cj> deter- 
mines the amount of energy transfer per unit time and 
unit volume. Using Eqs. (|A1|) and (|A2|1 one gets the 
following Langevine equation for the temperature T: 



dT 
~dt 













n + —(t> 


T 


r 





(A3) 



with coefficient r = c p p/-f. As stated before, £(t) de- 
scribes stochastic changes of temperature in time. Let 
us assume that these changes are such that their mean 
value is zero, 



(£(<)) = 0, 



(A4) 



whereas, for sufficiently fast changes, its correlator is 
equal to 



{Z(v)Z(t + At))=2D6(At). 



(A5) 



8 



The constants r and D define, respectively, the charac- 
teristic mean time for the temperature changes and their 
variance: 

(T(t)) = T + T(t = 0) cxp (-1} , (A6) 



(T 2 (t = oc)) 



-Dt. 



(A7) 



Thermodynamic equilibrium in such a situation means 
that for t >> t the influence of the initial condition 
Tit = 0) vanishes and the mean squared T has a value 
corresponding to the state of equilibrium. 

Eq. (|A3j> leads to the corresponding Fokker-Planck 
equation, 



df(T) 
dt 



d Id 2 
Q f K 1 f(T) + -—K 2 f(T), 



(A8) 



where the intensity coefficients are 

1 / r 
- T + 

T \ C p p 

n2 



D-2- 



K\{T) -- - j T„ 

K 2 (T) = 2DT 
Its solution is gamma distribution in 1/T = (3 



T, 



f(T) = 



r(a) \T 



where 



1 

7d' 



/' 



1 



T + — c 



(A9) 



(A10) 



(All) 



The mean value and the relative variance of (3 are, re- 
spectively (52j . 



W) = 

w 



T - —(, 
c p p 

1 

a 



tD. 



(A12) 
(A13) 



Temperature fluctuations in the form pre- 
sented here when applied to the exponential 
Boltzmann-Gibbs formula, Eq. |T]), lead to 

h q (E) = f °° exp(-E/T)f(T)d(l/T), i.e., to a Tsallis 
distribution, Eq. ([2]), with the nonextensivity parameter 
q equal to [l8( 



1 

Q 



q = l + - = l + rD 
and with the effective temperature 
T eff = T + — = T o + (g-l) 



r 
<i 



Depp ' 



(A14) 



(A15) 



In the case of no energy transfer, i.e., when = 0, one 
is left only with fluctuations and then T e ff — To 7 a s 
in [3]. Recently temperature fluctuations due to the 
volume fluctuations in statistical models were discussed 
as well and introduced in modelling of relativistic particle 
collision processes [53| . 



APPENDIX B: MULTIPLICITY DISTRIBUTIONS 
IN BOLTZMANN AND TSALLIS ENSEMBLES 

We shall now recapitulate some basic ideas concerning 
multiplicity distributions which result from, respectively, 
Boltzmann and Tsallis ensembles [2a|. 



1. Poisson multiplicity distribution 

This distribution arises in situations where, in some 
process, one has N independently produced secondaries 
with energies {-Ei,...,at} distributed according to a Boltz- 
mann distribution: 

g(Ei) = jexp(-^y, where A = (E). (Bl) 

The corresponding joint probability distribution is given 
by 



1 



9 ({-^1,..., at}) = — ]y exp . 



(B2) 



For independent energies {^i,„ jv} the sum E = 
^2^ =1 Ei is distributed according to a gamma distribu- 
tion, 

1 / E\ N ~ 1 ( E\ 
9n(E) = —777 m I T ) exp -- , (B3) 



X(N — 1)! V A 
the distribuant of which is 



N-l 



-{ (i-l)l VA 



EV 'expf—]. (B4) 



A 



Eq. (|B3p follows immediately either by using character- 
istic functions or by sequentially performing integration 
of the joint distribution (|B2|) and noticing that 



gN(E) = g N _x(E) 



E 



N - 1 



For energies such that 

N 



N+l 



J2 E i<E<J2Ei 



(B5) 



(B6) 



the corresponding multiplicity distribution has a Poisso- 
nian form (with (TV) = M-) 



P{N)=G N+ i(E)-G N (E) 



(N) 



N 



■exp(-(iV)). (B7) 



Therefore, whenever N variables {£ , i j ... i at i at+i} follow an 
exponential distribution (jBljl and satisfy condition (|B6| . 
then the corresponding multiplicity N has a Poissonian 
distribution (1571). 
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2. Negative Binomial multiplicity distribution 

This distribution arises when in some process N inde- 
pendent particles with energies {-Ei ... at} are distributed 
according to a Tsallis distribution, 



As before, for energies E satisfying the condition given 
by Eq. (|B6|) , the corresponding multiplicity distribution is 
equal to the well known Negative Binomial distribution 
(NBD): 



h({E 1 ,... tN }) = C 



N 



A 



- + 1-JV 



(B8) 

It means that there are some intrinsic (so far unspecified 
but summarily characterized by the parameter q) fluc- 
tuations present in the system under consideration. In 
this case we do not know the characteristic function for 
the Tsallis distribution. However, because we are dealing 
here only with variables {-Ei,...,jv} occurring in the form 

of the sum, E = > one can s ^ n l sequentially 

perform integrations of the joint probability distribution 
(IB8|) arriving at the formula corresponding to Eq. (|B3|1 : 



h N {E) = h N -x(E) 



E 



E 



N-l 



N-l (N-1)\X N 



N 



Y[l(i-2)q-(i-3)j 



l-(l-.)f 



(B9) 

- + 1-JV 



with distribuant given by 



N-l 

H N {E) = l-Y,H j {E) 



where 
H 3 {E) 



Ei- 



(j - 1)M 



jj[(i_2)g-(i-3)] 



i=i 



l-(l-«) 



E 



(BIO) 



(Bll) 

7+l-j 



P(N) = H N+1 (E) - H N (E) 



(B12) 



TV! 2-q 

N 



r ( 



i 



r(N + k) 



&-«>! 

m 



T(N + l)T(k) 



where 



N+k 



(B13) 



k 

Var(N) 



1 



E ~ 



l-{l-q) 



A 



(N) + (q - 1)(7V) 5 



(B14) 



In cases considered here q £ (1,2). In the limiting cases 
of q — > 1 one has fc — > oo and P{N) becomes a Poisson 
distribution whereas for q — ► 2 one has fc — > 1 and P(N) 
becomes the so called geometrical distribution. 
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